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The inverse matrix of some circulant matrices 

A. Carmona, A.M. Encinas, S. Gago, M.J. Jimenez, M. Mitjana 

Abstract. We present here necessary and sufficient conditions for the invert- 
ibility of circulant and symmetric matrices that depend on three parameters 
and moreover, we explicitly compute the inverse. The techniques we use are 
related with the solution of boundary value problems associated to second 
order linear difference equations. Consequently, we reduce the computational 
cost of the problem. In particular, we recover the inverses of some well known 
circulant matrices whose coefficients are arithmetic or geometric sequences, 
Horadam numbers among others. We also characterize when a general sym¬ 
metric circulant and tridiagonal matrix is invertible and in this case, we 
compute explicitly its inverse. 


1. Introduction 

Many problems in applied mathematics and science lead to the solution of linear systems having 
circulant coefficients related to the periodicity of the problems, as the ones that appear when using the finite 
difference method to approximate elliptic equations with periodic boundary conditions, see [I]. Circulant 
matrices have a wide range of application in signal processing, image processing, digital image disposal, 
linear forecast, error correcting code theory, see mm- In the last years, there have been several papers on 
circulant matrices that attend to give an effective expression for the determinant, the eigenvalues and the 
inverse of the matrix, see for instance 13 ®®- 

In this paper, we consider circulant matrices of type Circ(a, b,c ,..., c) and Circ(a, b, c,..., c,b). This 
type of matrices raise when dealing, for example, with finite differences for solving one dimensional elliptic 
equations, or when computing the Green function of some networks obtained by the addition of new vertices 
to a previously known one, see [7j. We give a necessary and sufficient condition for its invertibility. Moreover, 
as it is known, their inverse is a circulant matrix and we explicitly obtain a closed formula for the expression 
of the coefficients. 

For fixed n £ N*, we consider the vector space R" together with the standard inner product (•,•}. 
Moreover, we denote the components of the vector v S R” as Vj, j = 1,..., n, i.e., v = (iq,..., u„) T . 

As usual, 1 is the all ones vector and 0 is the all zeroes vector in R". In addition, e is the 1" vector 
whose first component is equal to 1 and 0 otherwise. 

The set of matrices of order n with real coefficients is denoted A'fn(R). Moreover, matrix i is the 
Identity matrix of order n and J is the matrix of order n whose entries arc all ones. 
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A matrix A = (a^) is named circulant with parameters a±,... ,a n if 


(1) 


or equivalently, 


ai 

a 2 • • 

a n 

a n 

ai • • 

* Oji—l 

a 2 

a 3 ■ ■ 

ai 


( 2 ) CLij ®l+(j-i)(modn) 

see 0[S]- 

Given a £ R", Circ(a) = Circ(ai,..., a n ) £ A4 n (R) is the circulant matrix with parameters a \,..., a n . 
Notice that Circ(e) = I and Circ(l) = J. 

Let r be a permutation of the set {1,... ,n} defined as, 


( 3 ) 


t(!) = f, r(j) = n + 2 - j, j = 2,...,n. 


We denote P r £ A4 n (R) the matrix with entries such that for any j = 1,... ,n, P T (j)j = 1 anc f 
Pij = 0, otherwise. It turns out that P T is invertible and it is satisfied Pf 1 = Pj = P T -i = P T . 

In a similar way we now define a r = P T a; that is, the vector whose components are (a T )\ = a± and 
( a T )j = a n + 2 -j , j = 2,... ,n. Thus, l r = 1 and (a r , 1) = (a, 1). Moreover, 

(4) Circ(a T ) = P T Circ(a)P T . 


For any a € R”, the matrices 


(5) Circ T (a) = P T Circ(a) = 


(6) Circ T (a) = Circ(a) P r = 


a-i 

a 2 

On 

a 2 

0-3 

Ol 


a 1 ■■■ 

On— 1 

ai 

On 

• • • a 2 

0"n 

On— 1 

• • • ai 

_ a 2 

ai 

• • • a 3 _ 


are named left-circulant and right-circulant with parameters a±,... ,a n , respectively. Both matrices are 
symmetric and for this reason matrices Circ T (a) were called circulant symmetric matrices in [3j. To avoid 
confusion, we will refer to this type of matrices with the notation introduced above. In addition, from 
identity (O we get Circ T (a) = Circ T (a r ) for any a € R". 

For the sake of completeness, we enumerate the properties of the circulant matrices that are relevant 
in our study. All the statements are easily shown to hold. 


Lemma 1.1. Given any a G R n , the following properties hold: 

(i) For any b € R" and a,/3 £ R, Circ(aa + /3b) = aCirc(a) + /3Circ(b). 

(ii) Circ(a) T = Circ(a r ). In particular, Circ(a) is symmetric iff a = a r . 

(iii) Circ T (a) = Circ(a) iff Circ T (a) = Circ(a). This equalities hold iff a = a r . 

(iv) Circ(a)l = (a, 1) 1. Moreover, if Circ(a) is invertible then (a, 1) ^ 0. 

(v) For any b £ R”, Circ(a)b = Circ(b T )a T and Circ(a)Circ(b) = Circ(b)Circ(a) = Circ(c r ), where 
c = Circ(a)b T = Circ(b)a. 

(vi) Circ(a) is invertible iff the linear system Circ(a)g = e is compatible. In that case, there exists a only 
solution g(a) that, in addition satisfies (g(a),l) = (a,l) _1 . Moreover, Circ(a) -1 = Circ(g(a)) and 
a r = a iff g(a) T = g(a). 
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(vii) Circ r (a) and Circ T (a) are invertible iff Circ(a) is invertible and, in that case is 
Circ T (a ) -1 = Circ T (g(a)) and Circ r (a ) _1 = Circ r (g(a)). 


One of the main problems in the field of circulant matrices is to determine invertibility conditions 
and, in this case, to calculate the inverse. From the previous Lemma, it turns out that for any a £ 1" the 
invertibility of matrices Circ T (a) and Circ T (a), as well as the computation of their inverses, can be deduced 
from the invertibility and the inverse matrix of Circ(a). Furthermore, the solution of this last problem can 
be reduced to study the compatibility of a certain linear system and the computation of its solution, when 
it exists. Moreover, the problem has been widely studied in the literature by using the primitive n-th root 
of unity and some polynomial associated with it, see 00 - Specifically, let u> = e « * be the primitive n-th 
root of unity. In addition, define for each j = 0,..., n — 1, the vector t ; = (l, ojf ... , up( n_1 )) £ R ra and 

n 

for any atM" the polynomial P 3 (x) = a j xJ_1 - Observe that to = 1 and for any a £ R n , P a (l) = (a, 1). 

j=i 

The following lemma provides a necessary and sufficient condition for the invertibility of Circ(a) and gives 
a formula for the inverse. 


Lemma 1.2. For any a £ R n , the following properties hold: 


(i) Circ(a)tj = P 3 (u) J )tj, for any j = 0,... ,n. In particular, detCirc(a) = J} P a (w fc ). 


k -o 


(ii) Circ(a) is invertible iff P 3 {u} 3 ) 7 ^ 0, j = 0,..., n — 1. In this case, Circ i (a) = Circ(h a ) where 

1 n —1 

( M 3 = - E oj-w-vp^)- 1 . 

n k = 0 


Let us notice that property (i) of the previous Lemma implies that all circulant matrices of order 
n have the same eigenvectors but different eigenvalues. Although the problem is completely solved, the 
computational complexity of formula (ii) for the determination of the inverse of a circulant invertible matrix 
grows with the order of the matrix, so it is not useful from the computational point of view. To illustrate 
this phenomena, we apply Lemma 1 1.2 1 to the cases n = 2,3, which can be directly solved without difficulty. 

For n = 2, uj = — 1 and for a given a £ R 2 , is P 3 (x ) = ai + a 2 X, and P a (—1) = ai — 02 - Thus, Circ(a) 
is invertible iff (ai + 02 ) (ai — 02 ) 7 ^ 0 and 


{h 3 \ 


( h *) 2 


1 / 1 


+ 


2 V ai + 02 ui 02 


2 V ai + a 2 ai — a 2 


ai 

of - a\ 
a 2 

of - a 2 


On the other hand, if n = 3, oj = |(—1 + iy/ 3), u> 2 = uj and for a given a £ R 3 , is P 3 (x) = ai + a 2 X + a^x 2 ] 
P a (w) = ai + 02 ta + a^UJ and P a (w 2 ) = P a (w). Thus, Circ(a) is invertible iff 

P a (l)|P a M ! 2 = (ai + o 2 + o 3 )|oi + o 2 w + o 3 w | 2 = al + a\ + o| - 3oia 2 a 3 7 ^ 0 


in this case, 




{h 3 \= 


{h 3 ) 3 — 


3P a (l)|P a (w )| 2 

1 

3P a (l)|P a (w )| 2 

1 

3P a (l)|P a (w )| 2 


(|P a M | 2 + 2P a (l)5R(P a ( W )) 
(\P 3 (oj)\ 2 + 2P 3 (m(ojP 3 (oj 
(|P a M | 2 + 2P a (l)5R(oJP a (a;; 


_Qj ~ Q2Q3_ 

of + a| + a| — 3aia2a 3 ’ 

a 2 — 0102 

af + of + a 3 — 3aia2a 3 ’ 
of — o 3 o 3 

af + of + a| — 3aia2a 3 
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We aim here to compute the inverse matrix of some circulant matrices of order n > 3 with three 
parameters at most. For instance, this kind of circulant matrices appears when computing the effective 
resistance and the Kirchhoff index of a network that comes from the addition of new nodes to a previously 
known one, see [7]. We significantly reduce the computational cost of applying Lemma [1.21 since the key 
point for finding the mentioned inverse matrix consists in solving a difference equation of order at most two. 


2. Matrices Circ(a, b, c,..., c) 

Given a, b, c £ M, let a(a, 6 , c) £ K n defined as a(a, b, c) = (a, b, c,..., c) T . Then, Circ(a, b, c,..., c) = 
Circ(a(a, b, c)). 

For any q £ M, we also consider the vector z (q) = (q n ~ 1 , q n ~ 2 , ■ ■ ■ ,q, 1) T £ R n . It immediately follows 

q n _ i 

that z T {q) = ( q n ~ 1 , 1, q ,..., q n ~ 2 ) T and (z(<?), 1} = -. We observe that the last identity also holds for 

q - 1 

q n — 1 

q = 1, since lim --= n = (z(l), 1). We remark that z(l) = 1 and z r (— 1 ) = z(— 1 ) = (—1,1,-1,..., 1) T , 

i <7 — 1 

for n even. 

Proposition 2.1. For any q £ R it is satisfied 

Circ(a(g, —1,0))z r (g) = [q n - l]e. 

Moreover, the following properties hold: 

(i) Circ(a(< 7 , —1, 0)) is invertible iff q n y^ 1, and the inverse matrix is 

Circ(a(< 7 , -1, 0 ))" 1 = (q n - l^Circ(z(g)). 


(ii) 


n 


The linear system Circ(a(l, — 1,0))h = v is compatible iff ^2 y i = 0. In this case, for any 7 G M the 

2—1 

unique solution of the linear system satisfying (h, 1 ) =7 is 


hj = — r 

n L 


7 ~ 12 ivi 

i= 1 


i=3 


3 = !,•••,«■ 


(iii) If n is even, the linear system Circ(a(—1, —1, 0))h = v is compatible iff (v, z(— 1)} = 0. In this case, 
all the solutions are given by 

n 

hj = (-1 y +1 + az(-l), j = 1 ,... ,n, a £ M 

*=i 


and hence 2 (h, 1 ) = (v, 1 ). 


The inverse matrix of a circulant matrix whose parameters are a geometric sequence can be computed 
as an application of the above result. 

Corollary 2.2. For any a,r £ M, the matrix Circ (ar 71 ^ 1 ,..., ar, a) is invertible iff a(r n — 1) 0. In this 

case, the inverse matrix is 

Circ(ar n_1 ,..., ar, a ) -1 = ( a(r n — 1)) 1 Circ(a(r, — 1 ,0)). 


The main result on the present section displays a necessary and sufficient condition for the invertibility 
of the matrices we are dealing with. Moreover, when the inverse exists, we can provide a simple and closed 
expression for its entries. 

Theorem 2.3. For any a,b,c £ R, the circulant matrix Circ(a, b, c,..., c) is invertible iff 

[a + b + (n - 2 )c] [(a - b) 2 + (l - (- 1 )") (c - b ) 2 ] y^ 0 
and, in that case, Circ(a, b, c,..., c) -1 = Circ(k(a, b , c)) where, if a ^ 2c — b 







The inverse matrix of some circulant matrices 


5 


j=l,. ■ ■ ,n, and 


( c—by 1 (a — c) n i c 

((a — c) n — (c — b ) n ) (a + b — 2 c) (a + b + (n — 2 )c) ’ 


fcj( 2 c — b, 6 , c) 


nine 


n — 1 
2 (c — 6 ). 


(J-1) 

n(c — 6 ) ’ 


J = 1, - - 


n. 


PROOF. Define a(n; a,b,c) = [a + 6 + (n — 2)c] [(a — b ) 2 + (l — (—1)”) (c — 6) 2 ] , thus a(n; a, 6, 6) = 0 iff 

a = b or a = (1 — n)b. Moreover, when c ^ 6, if n is odd, then a(n; a, 6, c) = 0 iff a = —6 — (n — 2)c, while if 

n is even, then a(n; a, 6, c) = 0 iff either a = b or a = —b — (n — 2)c. 

According to (iv) in Lemma [TTT1 (a(a, b, c), 1) = a + b + (n — 2)c ^ 0 is a necessary condition in order 
to ensure the invertibility of Circ(a(a, &, c)). So, from now on we will assume that this condition holds. 
Moreover, from Lemma [1.11 fvi). Circ(a(a, b, c)) is invertible iff the linear system Circ(a(a, b, c)) h = e is 
compatible, and in that case there is a unique solution that, in addition, satisfies (h, 1} = (a(a, b, c), 1) ■ . 

Therefore, since Circ(a(a, b, c)) = Circ(a(a — c, & — c, 0)) + cJ, if h £ R ra fulfills Circ(a(a, 6, c)) h = e, 
then Circ(a(a — c, b — c, 0))h = e — c(a(a, b, c), 1) -1 1. 

Reciprocally, if h £ K" satisfies Circ(a(a — c, b — c, 0))h = e — c(a(a, 6, c), 1) —1 1, then 

Circ(a(a, 6, c))h = Circ(a(a — c, b — c, 0))h + cJh = e + c[(h, 1) — (a(a, b, c), 1) _1 ] 1 


and hence, h is a solution of Circ(a(a, 6 , c)) h = e iff (h, 1) = (a(a, 6 , c), 1) 1 . 
Consequently, we have shown that 


Circ(a(a, b, c))h = e iff Circ(a(a — c, b — c, 0))h = e — c(a(a, 6, c), 1} 1 1 
and moreover (h, 1} = (a(a, b, c), 1) _1 . 

If c = b, then (a(a, b, b), 1} = a + (n — 1)6 and Circ(a(a — c, 6 — c, 0)) = (a — 6)1. Thus, the system 
(a — 6)lh = e — 6(a + (n — 1)6) is compatible iff a ^ 6 and then 


h = 


which implies (h, 1 ) = 


(a — 6 ) (a + (n — 1 ) 6 ) - 
(a + (n — 1 ) 6 ) — bn 1 


(a — 6 ) (a + (n — 1 ) 6 ) a + (n — 1)6 

1 6 


(a + (n — l) 6 )e — 61 , 

= (a(a, 6 , 6 ), l) -1 . Notice that, 


hi(a, 6 , 6 ) = 


hj(a 1 6 , 6 ) = — 


— 6 (a — 6 ) (a + (n — 1 ) 6 ) ' 

6 


(a — 6 ) (a + (n — 1 ) 6 ) 


j = 2,...,n. 


d — C 

For the case c ^ 6 , we consider q = --, then a (a — c, 6 — c, 0) = (c — 6 )a(g, —1,0) and, as a consequence, 

c — b 

the system Circ(a(a — c, 6 — c, 0)) h = e — c(a(a, 6 , c), 1) -1 1 is equivalent to 


Circ(a(< 7 , —1, 0))h 


-- -7 - 7 - 7 —r- ((a + b + (n — 2 )c)e — cl). 

(c- 6 )(a + 6 +(n- 2 )c) V v v ’’ ) 


If h is a solution of the previous system, then 


(a + 6 — 2 c) 

(c — 6 ) (a + 6 + (n — 2 )c) 


(Circ(a(g,-l,0))h,l) = (h, Circ(a T (g,-1,0))l) 
<a(g,-l, 0 ),l)<h,l)= (a ( c b _~ b) 2c) <M), 
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which implies that if a + b — 2c y^ 0, then (h, 1} =-= (a(a, 6 , c), 1) 1 . So that, if c b and 

a + b + (n — 2 )c 

a + b— 2c y^ 0; i.e. q 1, then Circ(a(a, b , c))h = e iff Circ(a(< 7 , — 1, 0))h = --r-r- ( (a + 

(c — 6 ) (a + b + (n — 2 )c) V 

b + (n — 2 )c)e — clj . 

For n odd or n even but b y^ a; i.e. g” y^ 1, according to Proposition 12.11 (i), Circ(a(< 7 , —1, 0)) is 
invertible and in addition, 


h = (c— b){q n — l)(a + 6 + (re — 2)c) Circ ^ z(< ^ > ) (< a + 4 + <» “ 2 >+ - + 

= (o- t )fa»-l)(I + fc + (n- 2 )c) ( ( ° + 6 + ( " “ 2)CM,) “ C(2(g) ' 1( 
If n is even and b = a, then <7 = —1, a + b + (n — 2)c = 2(6 — c) + nc and 
((26+ (n — 2 )c)e — cl,z(— 1 )) = — ( 2(6 — c) + nc) y^ 0 . 


1 . 


Consequently, the system 


Circ(a(—1, —1, 0))h 


-- —7 ---—x ((a + b + (n — 2 )c)e — cl) 

(c-6)(a + 6+(n-2)c) V v ) 


is incompatible, thus Circ(a(a, 6 , c)) has no inverse. 

If c y^ 6 and a + 6 — 2c = 0, i.e. <7 = 1, then a + 6 + (n — 2)c = nc and the system 


Circ(a(l, —1, 0))h = 1 (ne - 1) 

v 2 n(c — b) 

has solution. Moreover, by Proposition 12.II liil. if 


/ii 


n — 1 ' 

2(c — 6) J ’ 


h o = 


1 

Lnc 


n — 1 ' 
2 (c-b). 


(n + 1 - j) 
n(c — 6 ) 


j = 2, 


n 


then, h is the only solution of the linear system that satisfies (h, 1 ) = —. 

nc 

In all cases it suffices to take k = h T . □ 

Notice that for a = 2c — 6 in the previous theorem, Circ(a(2c — 6 , 6 , c) is invertible iff c(c — 6 ) y^ 0 
and then, the entries of vector h(2c — 6 , 6 , c) are the elements of an arithmetic sequence. Thus we can give 
a characterization of the inverse matrix of a circulant matrix whose parameters are in arithmetic sequence. 
Of course, our results coincide with those obtained in §5] and, for left-circulant matrices with parameters in 
arithmetic progression, in [3j. 

Corollary 2.4. For any a, b € R, the matrix Circ(a, a+ 6 ,..., a+(n— 1)6) is invertible iff (2a+(n— 1 ) 6)6 y ^0 
and in this case, 


Circfc 


In particular, for any m £ Z such that 2m + n / 1, the matrix Circ(m, m + 1,..., m + n — 1) is invertible 
and its inverse is 


Circ(m, m + 1,..., ro + n — 1) 1 = 


9 —X- 77 J - -Circ(a(l, -1, 0)). 

+ (/m + n-1) n 
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3. Matrices Circ(a, b,c,...,c,b) 

For any a, b, c £ R, let b(a, b, c) £ R” defined as b(a, b, c) = (a, &, c,..., c, b). Then, Circ(a, b,c,... ,c,b) = 
Circ(b(a, b, c)) and b r (a, b, c) = b(a, b, c), since matrix Circ(a, b,c,... ,c,b) is symmetric. Regarding the case 
b(a, b, b) = a (a, b, b ), matrix Circ(a, b,b,... ,b,b) has been analyzed in the previous section, so from now on 
we assume c ^ b. The case c = 0 has been analyzed in (TO] under the name of symmetric circulant tridiagonal 
matrix, assuming the condition |a| > 2\b\ > 0; that is, that Circ(b(a, 6,0)) is a strictly diagonally dominant 
matrix. 

Notice that Circ(b(2, — 1, 0)) is nothing but the so called combinatorial Laplacian of a n-cycle. More 
generally, for any q £ M, Circ(b(2g, —1,0)) is the matrix associated with the Schrodinger operator on the 
cycle with constant potential 2(g — 1) and hence its inverse is the Green’s function of a n-cycle; or equivalenty, 
it can be seen as the Green function associated with a path with periodic boundary conditions, see HU. 
Since the inversion of matrices of type Circ(b(2q r , —1,0)) involves the resolution of second order difference 
equations with constant coefficients, we enumerate some of their properties. 

A Chebyshev sequence is a sequence of polynomials {Q n (x)}ne z that satisfies the recurrence 

(7) Qn+i(%) = 2 xQ n (x) - Q n -i(x), for each n £ Z. 

Recurrence 0 shows that any Chebyshev sequence is uniquely determined by the choice of the corresponding 
zero and one order polynomials, Q o and Q\ respectively. In particular, the sequences { T n and 
{U n }n=-oo denote the first and second kind Chebyshev polynomials that are obtained when we choose 
Tq(x) = Uo(x) = 1, Ti(x) = x , Ui(x) = 2x. 

Next we describe some properties of the Chebyshev polynomials of first and second kind that will be 
useful in the present work. See |12j for proofs and more details. 

(i) For any Chebyshev sequence {Q n }n =-^ there exists a, f3 £ R such that Q n (x) = aU n -i(x) + 
PU n - 2 (x), for any n £ Z. 

(ii) T_ n (x) = T n (x ) and U- n (x) = —t/ ra _ 2 (x), for any n £ Z. In particular, U-\{x) = 0. 

(hi) T 2n+1 (0) = U 2n+ 1 (0) = 0, T 2n { 0) = U 2n { 0) = (-1)”, for any n £ Z. 

(iv) Given n £ N* then, T n (q) = 1 iff q = cos(^), j = 0,..., whereas U n (q) = 0 iff q = 

cos T ), j = 1,... ,n. In this case, U n -i(q) = (-1) J+1 and U n +i{q) = (-l)A 

(v) T n (l) = 1 and I7 n (l) = n + 1, whereas T n (— 1) = (—l) n and U n {— 1) = (—l)”(n + 1), for any n £ Z. 

(vi) T n (x) = xU n -i{x) - U n - 2 {x) y T' n {x) = nU n - i(x), for any n £ Z. 

n 

(vii) 2(x — 1) J2 Uj(x) = U n+ i{x) — U n (x) — 1, for any n £ N. 

j=o 

Chebyshev recurrence m encompasses all linear second order recurrences with constant coefficients, 
see na. so we can consider more general recurrences. Let {H n (r, s)}^T 0 , where r, s £ Z and s ^ 0, the 
Horadam numbers defined as the solution of the recurrence 

(8) H n +2 = fH n+ i + sH n , Hq = 0, H\ = 1. 

Notice that for any n £ N*, H n ( 1,1) = F n , the n-th Fibonacci number , H n ( 2,1) = P n: the n-th Pell 
number , H n ( 1,2) = J„, the n-th Jacobsthal number and H n ( 2,-1) = ?7 n _i(l) = n. 

The equivalence between any second order difference equation and Chebyshev equations leads to the 
following result, see [T3i . Theorem 3.1] and [21 Theorem 2.4]. 

Lemma 3.1. Given r, s £ Z and s ^ 0, we have the following results: 

(i) If s < 0, then H n (r , s) = (y/~Cs) n ~ 1 U n -i{ 2 ^_ s ), n £ N*. 

(ii) If s > 0, then H 2n (r,s ) = rs n ~ 1 U n - 1 (l + ^), n £ N*. 

In particular, for any n £ N*, F 2n = U n _i( |), J 2n = 2 n ~ 1 U n _i (|), and P 2n = 2C/„_ 1 (3). In addition, 
H 2n {r,r) = r n U n -\(l + §) when r > 0 and H n (r,r) = (■ K / r ^r) n ~ 1 U n -i(^-) for r < 0. 
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In addition, for any ijgRwe denote by u (g), v(g) and w(g) the vectors in R” whose components are 
Uj = Uj- 2 {q), Vj = Uj-i(q) and Wj = Uj- 2 {q) + U n -j(q), respectively. 


Lemma 3.2. For any q € M", the following properties hold: 


(i) w T (g) = w(g) and (w(g), 1 ) 


Tn(q) - 1 
g-1 


Moreover, w(l) = nl. 


(ii) w (q) = 0 iff q = cos(^), j = 1 ,..., r^l- In this case, (u(g),l) = (v(g), 1 } = 0 . 

(iii) When n is even, then w 2 j-i(0) = 0 and w 2 j{ 0) = (— l )- 7-1 [l — (—l)v ], j = 1,..., 

(iv) When n is odd, then w 2 j-i(0) = (— l) 2 r i+ - 3 , j = 1,..., and w 2 jff)) = (—l )- 7-1 
(iv) When n is odd, then Wj(— 1) = (—l) J_1 (n + 2 — 2 j), j = 1,... , n. 


n 

2 ' 



n— 1 
2 


PROOF. w(g) = 0 iff U n -j(q) = —Uj- 2 (q) for any j = 1,. .. ,n and this equality holds iff U n -\{q) = 0 and 
U n - 2 (q) = —1. Moreover, U n -i(q) = 0 iff q = cos (^f), k = 1,..., n — 1, thus U n - 2 (q) = (—l) fc+1 , leads to 

U n - 2 (q) = -1 iff k = 2 j. □ 

Remark: The quotient is well defined for q = 1, because T n ( 1) = 1, U n (l) = n + 1, and 

T'(q) = nU n -i(q), using l’Hopital’s rule, lim(w(g), 1) = riU n -\(l) = n 2 . Moreover, for q = 1, is w(l) = nl 

g-s-l 

thus, (w(l), 1} = n 2 . 


Proposition 3.3. For any q £ R, 

Circ(b( 2 < 7 , -1, 0))w(g) = 2[T„(g) - l]e. 

and the following holds: 

(i) Circ(b(2g, — 1,0)) is invertible iff q^ cos (^ti), j = 0,..., l" 1 ^! and, 

Circ(b(2g, —1, 0 )) _1 = ^ ^ ^ Circ(w(<?)). 

(ii) If q = 1, the linear system Circ(b(2g, —1, 0))h = v is compatible iff (v, 1) = 0 in this case, for any 
7 E R the only solution satisfying (h, 1 ) =7 is given by 

1 n 
^ - 1 - % 

hj = -— > \j - i\(n - \i -j\)vi, j = 1 ,... ,n. 

n 2n 

i =1 

(iii) If q = cos(-^), j = 1,..., [2^1], the linear system Circ(b(2g, —1, 0))h = v is compatible iff 
(h,u (q)) = (h, v(g)) = 0 . 


Proof. To prove (i), notice that w (q) is the first column of the Green function for the Schrodinger operator 
for a n-cycle, or equivalently for a (n+l)-path with periodic boundary conditions, see [111 Proposition 3.12]. 


To prove (ii), it suffices to see that G = (gij), where (n 2 — l — 6\i—j\(n— \i— j|)), i,j = l,...,n 

is the Green function of the Combinatorial Laplacian of the cycle, see for instance m- The third claim (iii), 
comes from (ii) of Lemma 13.21 that states w(g) = 0. In addition, in this case, U n -i(q) = 0, U n - 2 (q) = —1 
and U n (q) = 1. Besides, vectors u(g) and w(g) satisfy 


2qui -u 2 -u n = —l- U n - 2 (q) = 0, 

-ui - u n - 1 + 2qu n = -U n - 3 (q) + 2qU n - 2 (q ) = U n -i{q) = 0, 


2qv\ -v 2 -v n =2q-2q- U n -i(q) = 0, 

-vi - v n -i + 2qv n = -1 - U n - 2 (q) + 2g[/ n _i(g) = 0, 

thus, Circ(b(2g, — 1, 0))u(g) = Circ(b(2g, —1,0))v(g) =0. 


□ 
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Next, the main result in this section is proved. We give necessary and sufficient conditions for the 
existence of the inverse of matrix Circ(a, 6 , c,..., c, 6 ) and we explicitly obtain the coefficients of the inverse, 
when it exists. 


Theorem 3.4. For a,b,c £ R, the circulant matrix Circ(a, 6 , c,.. 

PrM 

(a + 2b + (n — 3)c) [a — c + 2(6 — c 

j=i 


., c, 6 ) is invertible iff 



and, in this case 

where if a ^ 3c — 26 
9j [a, 6 , c) 


Circ(a, 6 , c,..., c, 6 ) 1 = Circ(g(a, 6 , c)), 

U j - 2 {q) + Un-j{q) _c_ . 1 

2(c — 6 )[T rl (( 7 ) — 1] (a + 26 — 3c)(a + 26 + (n — 3)c) ' 


with q 


c— a 
2(6 - c) 


whereas 


9j (3c — 26,6, c) = 


1 


12 n(c — 6 ) 


(n 2 - 1 - 6 (j - l)(n + 1 - j)) + 


1 

n 2 c’ 


j = l,...,n. 


PROOF. From claim (iv) of Lemma o a necessary condition for the invertibility of Circ(b(a, 6 , c)) is 
(b(a, 6 , c), 1) = a + 26 + (n — 3)c ^ 0, so, we will assume that this condition holds. Moreover, claim (vi) 
in the same Lemma states that a necessary and sufficient condition to get Circ(b(a, 6 , c)) invertible is the 
compatibility of the linear system Circ(b(a, 6 , c))g = e, and in that case there is an only solution that satisfies 
(g, 1) = (b(a, 6 , c), l) -1 . 

As in Theorem 12.31 

Circ(b(a, 6 , c))g = e iff Circ(b(a — c, 6 — c, 0))g = e — c(b(a, 6 , c), 1) _1 1 

and moreover, (g, 1 } = (b(a, 6 , c), l) -1 . 

Since b(a — c, 6 — c, 0) = (c — 6)b(2g, —1,0), the linear system 

Circ(b(a — c, 6 — c, 0))g = e — c(b(a, 6 , c), 1) -1 1 

is equivalent to system 

Circ(b(2,, -1,0))g = (c _ t)(a + a + („- 3 )c ) ((“ + *+("- ~ <*) ■ 


If g is a solution of the above system, then 
(a + 26 — 3c) 


(c — 6 ) (a + 26 + (n — 3)c) 


= (Circ(b( 2 < 7 , —1,0))g, 1) = (g, Circ(b(2g,-1,0))l) 
= (b( 2 g, - 1 , 0 ), l)(g, 1 ) = (a + 26 ~ 3c) (g, 1 ). 


As a consequence, if a+26—3c ^ 0 then (g, 1 ) = 


1 


(c- 6 ) 

= (b(a, 6 , c), 1) . Under this assumption; 


a + 26 + (n — 3)c 

that is, if a ^ 3c — 26 or equivalently q 1, then Circ(b(a, 6 , c))g = e iff 

Ci,c(b(2„, -1,0))g = + 26 + (n - 3)c)e - cl). 
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M r 


— c + 2(6 — c) cos 


(¥) 


In addition, if 

j = 1 

Using claim (i) in Proposition 13.31 Circ(b(2g, — 1,0)) is invertible, and 

1 


P 0, then q p cos , for any j = ■ 


g : 


2(c - 6) (a + 26 + (n - 3)c) [T n (q) - 1] 
1 

2(c - 6) (a + 26 + (n - 3)c) [Tffq) - 1] 


Circ(w(g)) ^(a + 26 + (n — 3)c)e — cl^ 
^(a + 26 + (n — 3)c)w(g) — c{w(q), 1) 1^ . 


If there exists j = 1,..., such that a — c + 2(6 — c) cos = 0, i.e. q = cos , then, 

statement (ii) in Lemma 13.21 ensures 

((a + 26 + (n — 3)c)e — cl, v(g)) = (a + 26 + (n — 3)c)ui(g) = a + 26 + (n — 3)c p 0 

so, by claim (iii) in Proposition l3.3l the linear system Circ(b(a, 6, c))g = e is incompatible and, Circ(b(a, 6, c)) 
is not invertible. 

When a = 3c — 26, this is q = 1, then a + 26 + (n — 3)c = nc and system 

Circ(b(2, -1, 0))g = ^ (ne - 1) 

is compatible. Moreover, using claim (ii) in Proposition 13.31 the vector g € K" whose components are given 
for any j = 1,..., n by 

11 I 11 

S1 = ^ 0 - U" - {i ~ 1)1 + - b) g lj 

is the only solution of the system satisfying (g, 1) = —. Last, we only have to take into account that 

nc 

n n n 

E 1 3 - i\( n - I* - j\) = 77 - 1 ), for any j = 1 ,... ,n. □ 

i =1 6 


The case a = 3c — 26 in the above theorem, involves the Green function of a cycle. Cases related to 
this, raise as application in the analysis of problems associated with this combinatorial structures. 


Corollary 3.5. For a given a, 6 € R, matrix 

A = Circ(a, a + b(n — 1), a + 26(n — 2),..., a + jb(n — j),..., a + b{n — 1)) 
is invertible iff (6a + 6(n 2 — 1))6 p 0 and, 

A ' 1 = »y6»+tV — iP -^ 0 ^ 2 '- 1 '” 1 )- 

Corollary 3.6. For a given o, 6 £ K, the following results hold: 

(i) If n = 1 mod(4), then A = Circ(a, a, 6, 6, a, a,..., a, a, 6, 6, a) invertible iff (a — b)(a{n + 1) + 6(n — 
1)) p 0 and then 

A- 1 = —!_Circ(b(0,1,0))- . 2 ( a + b ) - - J 

a — 6 (a — 6) (a(n + 1) + 6(n — 1)) 

(ii) //n = 2mod(4), then A = Circ(ptpa, 6, ..., ptp 6, a) is invertible iff {a — b){a{n-\- 

1) + 6(n — 1)) p 0 and then 

A' 1 = -Circ(b(0,1, 0)) - -- -TV J 

a — 6 (a — 6) (a(n + 1) + b(n — 1)) 

(iii) If n = 3 mod(4), f/iea A = Circ(6, a, a, 6, 6,..., a, a, 6, 6, a, a) is invertible iff (a — 6) (a(n + 1) + 6(n — 
1)) p 0 and then 
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A " 1 = 


1 


-Circ(b(0,1, 0)) — 


2 (a + 6 ) 


a — b v v ’ ’ " (a — b) ( a(n + 1 ) + b(n — 1 )) 

(iv) When n is odd, then 

A = Circ(a + nb, a — (n — 2 )b,... , a + (—1 y~ l {n + 2 — 2 j)b,... , a — (n — 2)6) 
is invertible iff b(an + b) 0 and then 

We end up this paper by deriving the inverse of a general symmetric circulant tridiagonal matrix, 
without assuming the hypothesis of diagonally dominance. Notice the difference between our result and the 
methodology given in [ 10 ] • 

Corollary 3.7. For a, b € R, b 0, the circulant matrix Circ(a, b, 0,..., 0, b) is invertible iff 

rVl - '27rr 


fl [a + 26cos(//)] 


0 


j =o 


and, in this case 
where 


9 j(a,b, 0 ) = 


Circ(a, b, 0,..., 0, b) 1 = Circ(g(a, b, 0)), 

(_1H Mi) + 


26[1 - (~l) n T n (f b )\ L 


j = l,...,n. 


Notice that the diagonally dominant hypothesis \a\ > 2|6| clearly implies that a + 26 cos 0 for 

any j = 0 ,..., n. 
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